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Abstract 

The spin-one-half Falicov-Kimball model is solved exactly on an infinite- 
coordination-number Bethe lattice in the thermodynamic limit. This model 
is a paradigm for a charge-transfer metal-insulator transition where the occu- 
pancy of localized and delocalized electronic orbitals rapidly changes at the 
metal-insulator transition (rather than the character of the electronic states 
changing from insulating to metallic as in a Mott-Hubbard transition). The 
exact solution displays both continuous and discontinuous (first-order) tran- 
sitions. 
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I. INTRODUCTION 

The Falicov-Kimball Model (FKM)§ is the simplest model for a charge-transfer metal- 
insulator transition. In this model there are two types of electronic states: (i) localized 
f- or d-orbitals which have a negligible overlap with neighboring electronic orbitals but 
possess a strong on-site Coulomb repulsion (usually taken to be infinite), and (ii) delocalized 
conduction-band orbitals in which the Coulomb repulsion between two conduction electrons 
is neglected. The only "dynamic" Coulomb interaction that is included is the Coulomb 
repulsion between a conduction electron and a localized electron that occupy the same 
lattice site. The metal-insulator transition takes place when there is a rapid change in the 
thermodynamic occupation of the electron levels as a function of temperature (or some 
other thermodynamic variable such as pressure). Hence, an insulator (or semiconductor) 
which has most electrons lying in the localized states rapidly changes it's character to a 
metal as the electronic charge transfers from the localized levels to the conduction band. 
This transition is similar to the liquid-gas phase transition, in which the density changes 
abruptly at the first-order transition temperature or smoothly varies when the temperature 
is above the critical point. Here it is the density of electrons in the conduction band that 
changes at the transition. 

It was originally thought that this model described the most important physics behind 
the metal-insulator transitions in a variety of transition-metal and rare-earth compounds. 
These materials (such as SmB 6 , V2O3, Ti 2 03, NiS, etc.) exhibit a variety of metal-insulator 
transitions with either continuous or discontinuous changes in the conductivity occurring as 
the temperature is varied. This opinion, however, was not shared by everyone, since most of 
these materials simultaneously exhibit structural and metal-insulator phase transitions, and 
because the FKM ignores all effects arising from the hybridization between the localized and 
the conduction-band orbitals (such as the Kondo effect, and the screening of the localized 
moments). It remained unclear whether it was the electronic system, modeled by the FKM, 
that was driving the transition, or whether it was driven by other (electron-phonon or 
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hybridization) effects. 

Furthermore, there are competing scenarios for metal-insulator transitions. In the Mott- 
Hubbard scenario,! it is the electron-electron correlations that change the character of the 
electrons (within a single band) from an insulator to a conductor. Whereas in the Anderson 
picture! it is disorder that produces states that are localized and delocalized (within a single 
band), and the metal- insulator transition takes place by adjusting the Fermi level between 
the localized and delocalized states. 

Recently, however, a series of experiments! have been performed on N1I2, which appears 
to be a paradigm for a charge-transfer metal-insulator transition as described by the FKM. 
In this material, the Ni ions donate one electron to each of the neighboring I ions, filling 
the Iodine p shell. The Fermi level lies within the d-band of the Ni ions. This material 
is an insulator because of the strong electron correlations within the d-band (which also 
lead to antiferromagnetic order at low temperatures). As the pressure is increased, the Neel 
temperature increases, and then disappears at the same point where the conductivity has a 
discontinuous increase by two orders of magnitude. Detailed x-ray diffraction studies showed 
that there was no structural phase transition occurring at the metal-insulator transition. 
Instead the experimental evidence points toward a transfer of charge from the I ions to 
the Ni ions, which quench the local moments (changing Ni ++ to Ni + ), and leave behind 
conduction holes in the Iodine p-band. Such a scenario is described by a FKM.! 

The theoretical side of the FKM has also been controversial. The original solution of 
the mo del! involved just a mean-field analysis of the interband Coulomb interaction. That 
approximate solution found both continuous and discontinuous metal-insulator transitions. 
Later, an approximate version of the Coherent Potential Approximation (CPA) was applied 
to the model by different groups and produced conflicting results:!^ the first-order metal- 
insulator transitions were obtained in one solution;! while only continuous transitions were 
obtained in the other solution! Hence it is important to be able to solve the FKM in an 
approximation- free fashion (and in the thermodynamic limit), to resolve the question of 
whether or not it contains discontinuous metal-insulator transitions. 



We present an exact solution of the simplest spin-one-half Falicov-Kimball model on an 
infinite-coordination Bethe lattice (in the thermodynamic limit). We find that the solution 
possesses both continuous and discontinuous metal-insulator transitions. In Section II, we 
describe the formalism used in solving the problem, and describe the numerical techniques 
employed. In Section III, we present our results, and end with our conclusions in Section 
IV. 

II. FORMALISM 

The FKM consists of two types of electrons: a localized (dispersionless) valence band 
and a conduction band (separated by an energy gap A). The Fermi level lies within the 
energy gap at zero temperature, so that all of the electrons a priori occupy the localized va- 
lence band, and the system is an insulator. Hence, the FKM was originally described within 
an electron-hole picture where one considers holes within the valence band. The valence 
holes have a direct Coulomb repulsion Ujj, which disfavors two holes occupying the same 
localized orbital. In addition, there is an interband on-site electron-hole Coulomb interac- 
tion {—U < 0) which drives the charge-transfer metal-insulator transitions. The Coulomb 
interaction is attractive, because the electron and hole have opposite charge. The resulting 
Hamiltonian is 

H = ]T[e(k) + A + y]<<A. 4AJlk' + C ff y,frfrf;..l). , (1) 

k,cr i : &(T f i 

where d\. a (d^ a ) is the creation (annihilation) operator for a conduction electron of wave 
vector k and spin a, e(k) is the dispersion of the conduction band (with bandwidth W), 
and f} a (fi a ) is the creation (annihilation) operator for a localized hole at lattice site i (the 
on-site energy of the localized hole is chosen as the origin of the energy axis). Since Uff is 
largeB for most localized levels, we choose the limit where — > oo and restrict the number 
of the f- holes per site to nj < 1. Since the d-electrons in the conduction band originate 
by thermal excitation from the f-band, their number is constrained to rid = n* in order to 
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conserve the total number of electrons. 

We have found it to be more convenient to study the model in the particle picture rather 
than in the original hole picture given in Eq. ([!]). We employ a particle- hole transformation 
(U — ► fl an d flU ~ * 1 — flU) and represent the kinetic energy in the localized basis 
to transform the Hamiltonian into 

H = - ]T tijd\ a d ia + E f J2 flU - fJ. Yl( d \aU + flU) 

+ UJ2 dld la flu + u ff J2 flu flu , (2) 



where d\ a (dio) is the creation (annihilation) operator for a conduction-band electron of spin 
cr at site i, is the hopping matrix between lattice sites % and j [which yields the band 
structure e(k)], Ef = U — A — ^ is the localized electron site energy measured from the 
middle of the conduction band, and Uff is the on-site Coulomb repulsion between f-electrons. 
A chemical potential /i is introduced to satisfy the constraint + n/ = 1 as Uff — > oo. 
This Hamiltonian has also been used as a model for intermediate valence problems,i as a 
thermodynamic model for an annealed binary alloy,0 as a simplified Hubbard model,0 and 
as a model for metamagnetism and anomalous magnetic response^ (when a magnetic field 
is added). 

We choose to solve the FKM in the infinite-coordination-number limit following the 
method of exact solution developed by Brandt and MielschS and expanded by Freericks and 
Zlatic.0 In the infinite-coordination-number limit, the local approximation becomes exact, 
implying that one can neglect the momentum dependence of the irreducible self energy 
and the irreducible vertex functions, but one needs to determine the frequency dependence. 
Hence, the problem on the lattice can be mapped onto a problem on a single-site, but coupled 
to an effective medium, which represents the dynamical information of all of the other sites 
of the lattice. The effective medium needs to be determined self-consistently in order to 
solve the lattice problem exactly. 

We begin with the local Green's function of the conduction electrons for each spin, 
G{iuj n ) = G n evaluated at the Fermionic Matsubara frequencies u n = (2n + l)7rfcsT, and 



express it explicitly in terms of the "bare" Green's function Go(iu n ) (which contains all of 
the dynamical information of the other sites in the lattice): 

_ l-n f n f 

where 



1 + gexp {(E f - fi - U)/k B T} JJ{1 - UG (iu n )} 



- -i 

-i 



(4) 



Here q = (2J + 1)/ (2 J + 1) is the ratio of the spin degeneracies without (Jo) and with (J) 
an f-electron at a lattice site, and we take q — |. (This is a different value from the original 
FKM work which took q = 2. We do not expect the results to depend strongly on the value 
of q.) The "bare" Green's function also satisfies the usual Dyson equation: 

GZ 1 (iu n ) = G- 1 -V n , (5) 

which can be viewed as a definition of the self energy Yi{iuj n ) = S n . The loop for determining 
the Green's functions is completed by evaluating the self-consistency equation by summing 
the momentum-dependent Green's function over all momentum [i.e., integrating over the 
noninteracting density of states D(e)], 

f°° Die) 

G n = de - ■ — — , (6) 

iu n + n - e(k) - h n 



yielding an exact solution of the model. 

The algorithm for solving the FKM is the same as that used for numerically solving the 
Hubbard model:ll! (i) Begin with an initial self energy E n (init.) (we chose either £ n (init.) = 
or we interpolated from a higher-temperature run); (ii) use Eq. (|6|) to find G n ; (iii) then use 
Eq. (^) to find Go(iu n ). (iv) Next determine nf from Eq. ([|); and (v) determine the new 
G n from Eq. (|3]). (vi) Finally use Eq. (|5]) with the new Green's function and the old "bare" 
Green's function to extract the new self energy and (vii) go back to step (ii) to repeat the 
iteration until convergence is reached. 
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Since we are interested in dynamical properties, we also need to solve for the retarded 
Green's functions on the real axis. We do this by first performing an imaginary-axis calcu- 
lation to find the filling of the f-electrons (n/), and then solving the analytically continued 
equations ([|), (|5|), and @ where the Matsubara frequencies are simply replaced by the real 
frequencies {iuj n — > uo + i5). Convergence of these equations is rapid under iteration on both 
the real and imaginary axes (less than one hundred iterations on average for a convergence 
of the self energy to one part in 10 9 on the imaginary axis and to one part in 10 3 on the real 
axis) . 

Note that the self-consistent equations for fixed rif, i.e., Eqs. (H), (1), and are 
identical to those employed in the CPA. However, nf is explicitly determined by Eq. @ in the 
infinite-coordination-number limit, while, in the CPA, it is determined thermodynamically 
by minimizing the following trial free energy,!'! FcpA[ n f\i as a function of nj\ 



/oo 
deA(e) ln(l + e-^ kBT ) 
-oo 



+ ( M _tf + A+y)(l-n / ) 



+ k B T 



1 — n t 

nfinrij + (1 — nf) In 



(7) 



Q 

which, for the conduction electrons, has the form of a noninteracting electron system with 
the noninteracting density of states -D(e) replaced by the interacting one, A(e) = — ^lmG(e). 
Here e is measured from the chemical potential [i. Since the conduction electrons are effec- 
tively noninteracting for fixed value of (i.e., the interacting density of states for fixed 
nf does not vary with temperature), the trial free energy in Eq. (|7]) is also exact. Naturally, 
the rif obtained in the CPA by minimizing the trial free energy agrees with the nf found 
within the Brandt-Mielsch formalism using Eq. (|J). 

We solve the model on an infinite-coordination-number Bethe lattice. As the coordination 
number Z increases, the hopping integral t between the nearest neighbor sites is scaled as 
t — > ^ in order to have a nontrivial kinetic energy, and the density of states for the 
noninteracting system becomes Wigner's semicircle with W = 4t*: 
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D(e) = 2^W^ 2 -e 2 - (8) 

We take t* as our energy unit (t* = 1). The integral that defines the local Green's function 
[Eq. (||)] can now be performed analytically, to yield 



G(z) = — 2 - sgn[lm(z)] V (9) 



with z = z + \i — £(2) and z an arbitrary complex number. 

In addition, one can also calculate susceptibilities for charge-density-wave order or spin- 
density-wave order. When these susceptibilities become infinite the system has a second- 
order phase transition to an ordered state. We will not describe in detail how to calculate 
such susceptibilities here (a discussion has appeared already0'0) , because the system never 
underwent any second-order phase transitions for all of the parameters we considered in this 
work. 



III. RESULTS 

We begin our discussion in the low-conduction-band-density limit {rid 0) which leads 
to the "excitonic phase" considered by Ramirez-Falicov-Kimball.B The excitons (bound 
electron-hole pairs which do not contribute to the dc-conductivity) can be shown to form 
when U is larger than a critical value U c . The critical value of U was determined by solving 
the single-exciton problem,S where the Hamiltonian in Eq. ([I]) with one hole and one elec- 
tron produces a bound state just below the conduction band. This occurs when the Green's 
function, evaluated on the real axis at the lower band edge, is equal to —1/U, or 

1 r 2 Die) 
U c J-2 e + 2 + id 

for the infinite-coordination- number Bethe lattice after using Eq. (|9|). This is the same 
critical value of U at which the interacting density of states splits into two bands as — ► 0.0 
However, it is not clear whether this simple criterion for exciton formation (based on the 



S 



= | - {l ± ^l + 2(l-2n d )UG(u) + U*Gi(uj)} , (11) 



single-exciton problem) is sufficient to create an excitonic insulator when the conduction- 
electron density is small, but finite. Thus we analytically determine the conduction-band 
Green's functions at finite temperatures and finite, but small electron density. We find the 
solutions separate into two regimes: a weak-coupling regime with U < 1 and a strong- 
coupling regime with U > 1. 

The functional form of the self-energy is expressed as a function of the local Green's 
function and the <i-electron concentration asll 
U 

~~ 2G(c 

where the sign is chosen for each frequency in such a way that the self-energy is analytic. 
Since we are interested in the insulating phase {n d — ► 0), we expand Eq. ( p|) about n d = 
by factoring out (1 + UG) from the square root. If U < 1 (single-band regime), the factor 
1 + UG{uj) never vanishes [from Eq. fll0|) 1, so we must choose the negative sign, in order to 
have a vanishing self energy as U — > 0. If U > 1 (two-band regime), we must choose the 
positive sign for u < —2 in order to ensure analyticity (i.e., the right sign for Im E), since 
1 + UG(uj) is negative there. Consequently, we choose the minus sign for U < 1 (single 
band regime) and the plus sign for U > 1 (two-band regime). Hence, in the limit where 
< 1, we have 



T,(u;n d -> 0) 



1+UG(lj) 

TT UdU -f TT / 1 

(12) 

[ 1 + l7G(w) ' ' 

Combining this with Eq. (^), we obtain the functional form of the local Green's function 
for each case. The resulting form is 

G non {. - U) n d UG non ^-U) itu<l 

[1 + UG nm {u - u)\^{u-uy-± 

1 1 

-— + n d G non (u + — ; t* -> v/r^t*) , if £/ > 1 , 

where we used the noninteracting form of the Green's function, found in Eq. (Q) with z - 
uj + id. Note that the Green's function for U > 1 represents only the split-off lower band. 
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G(u;n d -> 0) = < 



(13) 



When U < 1, as expected, the Green's function has just a perturbed form from the 
noninteracting Green's function which is shifted by U, G non (uj — U), and the gap to particle- 
hole excitations at T = remains as A. Thus the number of the conduction electrons rid 
at finite temperatures is activated as in a semiconductor with a fixed energy gap. On the 
other hand, when U > 1, the Green's function for the lower band becomes a band-narrowed 
(t* — > y/ridt*) noninteracting Green's function^ (plus a constant shift by — -j) and has a 
weight rid (the upper band has weight 1 — rid). The unit-charge condition (rif + rid = 1) 
implies that this lower band is always half-filled, and the energy of the system (in the limit 
T — > and rid fixed) becomes 



Einsulator ~ > (1 ~ Tld)Ef + 2 de 6 



—ImG(u;n d -»• 0, U > 1) 

7T 



= 17- A-2 + n d (A + 2-C/-^-^V^) , (14) 

where the factor of 2 in the integral arises from the spin degeneracy of the conduction band. 
We will see below that this strong-coupling phase is not an excitonic insulator as believed 
by Ramirez-Falicov- KimballB for any finite temperature. It is interesting to note, though, 
that this insulating limit [where all of the electrons lie in the localized states with an empty 
conduction band (rid = 0)] is noninteracting at T = 0, because there are no conduction 
electrons, and hence the localized electrons do not feel any Coulomb repulsion. The ground- 
state energy is Ei nsu i ator = U — A — 2. The gap to particle-hole excitations (at T = for the 
insulating phase) is A for U < 1, but decreases toward zero as A + 2 — U — ^ for f/ > 1. 
At the critical value of U where the gap closes, one can see from the insulator energy in 
Eq. ([14]) that the insulating phase is unstable if A + 2 — U — jj < (or equivalently if 



U > 1 + + J A(l + j)) since the ground-state energy is lowered for small, but nonzero 
n d - 

There is another phase of the FKM that is also noninteracting. It is the metallic phase 
(for large U) where the electrons (one per site) in the valence band are all promoted to the 
conduction band (rid — 1)- In this case, the conduction band feels no Coulomb repulsion, 
because there are no /-electrons to scatter them. Thus the system is characterized by the 
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half-filled noninteracting conduction band. The energy of this metallic phase is then 

E metal = 2 fdeeD{e) = -^-. (15) 

J-2 07T 

Therefore, at T = 0, there is a transition from an insulating ground state to a metallic 
ground state when E metai < E insu i ator or 

U>A + 2--^. (16) 

Surprisingly, there is a small region ofU, A + 2- J F >C/>l + f + y/A(l + f) and a small 
region of the gap energy < A < ^(1 — ^) 2 « 0.026923, where the ground state is neither 
metallic nor insulating. In this nontrivial region, there must exist either an intermediate- 
valence state, or a charge-density-wave-ordered insulator. Detailed studies in this regime 
will appear in a future publication. 

In this contribution, we are interested in examining the discontinuous phase transitions 
between states that are connected either to the insulating phase or to the metallic phase as 
T —>■ 0. So we choose the bare gap to be large enough (A = 1), in order to be sufficiently 
far from any intermediate valence or charge- density- wave ordered phases. We vary U and 
for each value of U calculate the thermodynamic properties of the system. We expect 
interesting behavior to occur for U close to the metal-insulator transition point at T = 0, or 
U ~ 2.15117 [from Eq. fll6|)1. We also expect simple semiconducting behavior (with a gap 
A) to occur for U < 1. 

In Fig. |l|, we present our numerical results of rid (vertical cLXlS^ clS db function of 1/T (hor- 
izontal axis) for different values of U, when A = 1. The conduction-band filling rid is plotted 
on a logarithmic scale so that the linear behavior in 1/T indicates activated carriers as in 
a semiconductor (see U = 0.5). On the other hand, the metallic limit obtained in Eq. ( |T6|) 
for U > 2.15117 also agrees with the numerical results, since the system remains metallic 
for all temperatures. Fig. [I] shows both discontinuous (U = 2.150, 2.140) and continuous 
(U = 2.120, 2.000, 1.700) metal-insulator transitions. In the discontinuous transitions, the 
conduction-electron concentration rid follows the metallic solutions at high temperatures but 
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drops to the insulating solution at a sharp transition temperature. (For U = 2.150, rid drops 
by two orders of magnitude as the temperature changes by ~ 10~ 3 .) This electron concen- 
tration can be viewed as an approximation to the electrical conductivity, if we assume that 
in a real material there is also static disorder (from defects, impurities, etc.) so that at low 
temperatures the relaxation time r approaches a constant and the conductivity is propor- 
tional to ridT. We also calculate the intrinsic conductivity of the FKM below, assuming that 
all of the scattering of the conduction electrons arises from the localized /-electrons. 

The conduction-electron density of states, A(uS) = — -ImG(ou) (where uj is measured 
from the chemical potential \x) provides additional information about the metal-insulator 
transitions. We calculate it by solving for the Green's functions on the real-axis. Fig. |2| 
plots A(u) for some of the representative cases from Fig. [I] at various temperatures: (a) 
the metallic regime (U = 2.160) where the half-filled lower band increases in size as the 
temperature decreases; (b) the discontinuous metal-insulator-transition regime where the 
corresponding density of states in the lower band discontinuously collapses to the insulating 
phase at a critical temperature T c (0.060 < T c < 0.065 for U = 2.150); (c) the continuous 
metal-insulator-transition regime where the lower band is continuously reduced as the tem- 
perature decreases; (d) and the semiconducting regime where n d has an activated behavior, 
and A{uj) displays a finite pseudogap from the perturbed single conduction band. Note that 
all of the temperature dependence of the interacting density of states arises from the temper- 
ature dependence of rid, since the FKM with fixed n d and n/ has a temperature-independent 
density of states.0 

The exact solution of the FKM also allows the optical conductivity to be calculated from 
the following formula:0 

<r(u) = a f°°du/ f°°deD(e) A(e, J) A(e, J + u) - + ^ , (17) 

where f(oo) is Fermi distribution function and <To gives the conductivity unit. Note that this 
original formula was derived for the hypercubic lattice. It relies on two assumptions: the first 
is that the vertex corrections for the conductivity vanish in the infinite-dimensional limit, 
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and the second is that the derivative of the electronic band-structure (with respect to the 
crystal momentum) can be calculated. The former holds on a Bethe lattice, but there is no 
well-defined crystal momentum, and hence no well-defined bandstructure on a Bethe lattice. 
We believe (but have not been able to show explicitly) that the conductivity for the Bethe 
lattice should have the same form as on the hypercubic lattice except for some constants 
of order unity (which are absorbed into the definition of <7o). We use this assumption in 
calculating the dc conductivity <Jdc which is then found from Eq. ([17]) in the limit u — > 0. 

Substituting the spectral function into Eq. ( |17|) and assuming the self energy has neg- 
ligible frequency dependence near the Fermi level produces the following limiting form for 
the dc conductivity with T = and rid —> (i.e., for the insulating phase): 



a dc (T = 0;n d ^0) = , (18) 

which is proportional to the product of the intrinsic relaxation time (~ 1/ImS) and the 
density of states at the Fermi surface. Calculating the self energy E(0) from Eqs. flT2|) and 
(|T3|) with the proper chemical potential for each case (fi = U — 2 + 5fi for U < 1, and 
H = —jj for U > 1), we obtain 

( «°(i-uy ifc/<1 

(19) 



o~dc{T = 0;n d ^ 0) 



if U > 1 . 



I 4tt 2 (C/-1)V^ 

Thus the intrinsic Ode (which is obtained from the pure electronic system of the model) in 
the insulating phase at T = actually diverges as nd — > 0! Even in the strong-coupling 
regime (U > 1), the electron-hole excitations do not bind to form an excitonic insulator! 
(The dc conductivity does diverge more slowly (~ ^fjjj) than in the weak-coupling regime 
(~ as nd — > 0, though. ) Therefore, the excitonic phase considered by Ramirez-Falicov- 
KimballB does not exist on the infinite-coordination-number Bethe lattice. Moreover, the 
divergence of the intrinsic Od c occurs because the relaxation time r increases more rapidly 
than the density of states at the Fermi surface decreases when nd — > [see Eq. (|18|)1. But 

13 



in a real material the relaxation time can never diverge because there always exists some 
static disorder which forces the relaxation time to approach a constant at low temperature. 
Hence, the dc conductivity will approach zero as rid — > in any real material. In Fig. |3|, 
we present our numerical results of the intrinsic dc conductivity as a function of 1/T for 
the corresponding values of U in Fig. [!]. The conductivity in Fig. |3], with moderate rid 
appears to be proportional to the electron concentration shown in Fig. [I], but cr^c starts 
to diverge as temperature decreases and becomes small enough for the limiting form in 
Eq. © to hold. (For example, see the cases U = 2.000 and 1.700. For U = 0.500, a dc is 
always in the low-density limit. ) Thus, the relevant dc conductivity for a real material is 
approximated better by n^r with a constant relaxation rate r, rather than using the intrinsic 
dc conductivity! 

Finally, we examine the discontinuous metal-insulator transitions in more detail to show 
that they are indeed first-order phase transitions. We do this by employing the CPA for- 
malism to calculate the trial free energy as a function of n& (or equivalently nj = 1 — rid) 
at different temperatures. The minimum of the trial free energy determines the thermody- 
namic electron density. We first solve the real-axis self-consistent equations for each value 
of < nd < 1 to find the interacting density of states A(u). [Note that this A(u) at fixed 
nd is independent of temperature. 0] Then we evaluate the exact form of the free energy in 
Eq. (^) as a function of and repeat the process for different temperatures. We present 
our results near the critical temperature in Fig. |], where U = 2.150 and A = 1. The free 
energy has a double minimum near the critical temperature (0.060 < T c < 0.065) and as 
the temperature is lowered, the conduction electron density discontinuously changes as the 
global minimum switches between the two local minima (indicating a first-order transition). 
At T c , where the two minima are degenerate, the system exhibits phase coexistence between 
the insulating and metallic phases. In the region where the change in is continuous, the 
free energy does not have multiple minima, but rather the minimum of the free energy varies 
smoothly with nd as the temperature is changed. We also verified that the minimum of the 
CPA form of the free energy agrees with the form for the free energy determined by Brandt 
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and Mielsch.0 



IV. CONCLUSIONS 

In conclusion, we have exactly solved the spin-one-half Falicov-Kimball model on an 
infinite-coordination-number Bethe lattice, which is shown to have both continuous and 
discontinuous (first-order) charge-transfer metal-insulator transitions. By being able to solve 
the model exactly, we have clarified the theoretical controversy of the model and have proven 
that the model does display first-order metal-insulator transitions. The simplicity of the 
model, based on only the electronic system (which has both a localized and a conduction 
band), emphasizes the fact that the electronic system itself (the Coulomb interaction between 
a conduction electron and a localized electron) can cause dramatic discontinuous charge- 
transfer metal-insulator transitions without requiring other effects (such as phonons). We 
expect our results to continue to hold in three dimensions and to have applications to real 
materials such as NQ2. 

We also found that the intrinsic conductivity (determined by the scattering of the d- 
electrons off the f-electrons) actually diverges for the "insulating" phases with — > 
because the relaxation time grows faster than the density of states at the Fermi level de- 
creases. In a real material, the conductivity will go to zero as — > though, because the 
relaxation time is bounded by the scattering off of impurities. We also discovered a small 
region of parameter space that possesses intermediate-valence or charge- density- wave order. 
Further studies of this region are currently underway. 
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FIGURES 

FIG. 1. Number density of the conduction electrons rid, (on a logarithmic scale) plot- 
ted as a function of t*/T for various values of U, where A = l.Ot*. Four different 
regimes are shown: (a) the metallic regime for all T (U = 2.160**, 2.155**); (b) the 
discontinuous metal-insulator-transition regime (U = 2.150**, 2.140**); (c) the continuous 
metal-insulator-transition regime (U = 2.120**, 2.000**, 1.700**); and (d) the semiconducting 
regime {U = 0.500**). 

FIG. 2. Density of states A(uj) at different temperatures for some representative cases of Fig. |l]: 
(a) U = 2.160**; (b) U = 2.150t*; (c) U = 2.120**; and (d) U = 0.500** when A = 1.0**. Here the 
energy u is measured from the chemical potential fi (i.e., the Fermi level lies at u = 0). 

FIG. 3. The intrinsic dc-conductivity Ode of the Falicov-Kimball model as a function of t*/T 
for the corresponding values of U shown in Fig. ||. Note how the intrinsic conductivity diverges for 
low temperature and low electron concentration, as described in the text. 

FIG. 4. Free energy as a function of the number density of the conduction electrons near the 
first-order transition temperature. The parameters are A = 1.0** and U = 2.150**. The free 
energy has two local minima with the global minimum switching between these local minima as 
temperature varies through the critical temperature (0.060 < T c < 0.065), indicating a first-order 
transition. 



18 




19 



0.4 
0.3 
0.2 
0.1 
0.0 
0.3 
0.2 
0.1 
0.0 
0.3 
0.2 
0.1 
0.0 
0.3 
0.2 
0.1 



0.0 



I ' 

I (a)U=2.160t* 








(i 1 1 1 — 31 ^ 


\^ 

1 ^ 


_ (b)U=2.150t* 


~*~ 


Lh ^ — h 


\ \ " 

^ \ : 


1 (c) 11=2.1 20t* 


1 p — > — 






/ 

/ 


\ 
\ 

\ - 

\ ■ 

1 1 1 


(d) U=0.500t 

// 
i i I'l 


^^^^^ 

\ : 



T=0.070t 
T=0.065t* 
T=0.060t* 
T=0.055t* 



T=0.080t 

T=0.075t* 

T=0.070t* 

T=0.065t* 

T=0.060t* 

T=0.055t* 

T=0.400t* 
T=0.200t* 
T=0.150t* 
T=0.100t* 
T=0.070t* 



-2.0 



0.0 



2.0 



4.0 



CO 



20 




21 





1 






(a) T=0.070t* 


/ 


1 


^ i 

(b) T=0.065t* 


/ - 


1 


^ [ ~~ r 


- 


1 


^ 1 

i 





0.0 0.2 0.4 

n. 



0.6 



0.0005 

-0.0015 

-0.0035 

0.0010 

-0.0010 

0.0015 

-0.0005 
0.0035 
0.0015 
-0.0005 



22 



